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Abstract. - We generalize Einstein's master equation for random walk processes by considering 
that the probability for a particle at position r to make a jump of length j lattice sites, Pi{r) 
is a functional of the particle distribution function f{r,t). By multiscale expansion, we obtain a 
generalized advection-diffusion equation. We show that the power law Pj{r) oc ,f{r)"^^ (with a > 1) 
follows from the requirement that the generalized equation admits of scaling solutions (/ (r; t) = 
t~'' (j) {r /t'')). The solutions have a g-exponential form and are found to be in agreement with 
the results of Monte-Carlo simulations, so providing a microscopic basis validating the nonlinear 
diffusion equation. Although its hydrodynamic limit is equivalent to the phenomenological porous 
media equation, there are extra terms which, in general, cannot be neglected as evidenced by the 
Monte-Carlo computations. 



Introduction. — A standard procedure to describe the microscopic mechanism of a 
diffusion process, is to consider a test particle executing a random walk on some substrate. 
The idea goes back to Einstein who, in one of his celebrated 1905 articles [1], showed 
how the diffusion equation follows from a mean-field description written in terms of the 
probabilities that the particle performs elementary displacements at each time step. The 
distribution function /(r, t), that is the probability that, given the particle was initially at 
r = at t = 0, it will be at position r at time t (for t large compared to the duration of 
an elementary displacement) is obtained as the solution to the Fokker-Planck equation for 
diffusion, and one finds that, in the the long-time limit, f{r,t) is Gaussianly distributed in 
space [2]. This result had been known since Pick's law was established for diffusion; what 
was new in Einstein's work was the microscopic content, in particular the expression of the 
diffusion coefficient in terms of the particle velocity autocorrelation function, a form that 
was further generalized to the general class of transport coefficients known since the 1960's 
as the Green-Kubo coefficients [3]. 

The classical diffusion equation has been extensively used and successfully applied to 
a large class of phenomena (ranging from particle dispersion in suspensions to diffusion of 
innovations in social networks) and is indeed applicable as long as the system responds 
linearly to a change in the quantity that is being transported. But the linear response 
hypothesis doesn't hold in more complicated situations such as when there is an interactive 
process between the particle and the substrate, or in heterogeneous media. In the 1930's, 
a nonlinear diffusion equation was proposed on a purely phenomenological basis devised in 
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particular to describe diffusive transport in porous media, hence the name porous media 
equation [4] 

p{r;t)=D^nr;t), (1) 

where D is the diffusion coefficient. This equation, when generahzed with an advective term, 
has a g-Gaussian solution [5] and exhibits the interesting feature that the scaling (r^) cx H 
can be non-classical (7 7^ 1). ^ It was not until the f990's that a more fundamental 
basis was proposed for the (generalized) porous media equation using various statistical 
mechanical approaches. The reason for the various approaches can be found in the variety 
of problems where non-classical (non-Gaussian) distributions are observed: transport in 
porous media, viscous fingering, information diffusion in social networks or in the internet, 
financial market distributions, ... . The proposed approaches use formulations such as the 
generalized entropy [6], the Langevin equation [7,8], the master equation [9], the nonlinear 
response [10], the escort distribution [11], or the generalized generating hmction [12]. 

Generalized master equation. — Here we use the microscopic approach by going 
back to Einstein's original derivation based on the random walk. For simplicity consider a 
one-dimensional lattice where the particle hops to the nearest neighboring site (left or right) 
in one time step, a process described by the discrete equation 

n{r; t + 1) ^ (-n{r + l;t)+ £_+n{r - 1; t) , (2) 

where the Boolean variable n{r;t) = {0,1} denotes the occupation at time t of the site 
located at position r and ^± is a Boolean random variable controlling the particle jump 
between neighboring sites (^+ + i.^ = 1)- The mean field description follows by ensemble 
averaging Eq. (2). With {n{r;t)) = f{r;t) and (^j) = Pj (using statistical independence of ^ 
and n), and extending the possible jump steps over the whole lattice, one obtains Einstein's 
master equation [1] 

+00 

f{r;t + 6t) = ^ P, (r - j6r; t) f(r - jSr; t) , (3) 

j = -oo 

where Pj{i) denotes the probability that the walker at site £ make a jump of j sites ^. Using 
the normalization: 1 = Pq {£; t) + J2j^o ^3 ^) ' Ell- (3) takes the form of a Boltzmann like 
difference equation 

fir;t + St) - f (r; t) = ^[P, (r - jSr; t) f (r ~ jSr; t) - P, (r; t) f (r; t)] , (4) 

which simply describes the rate of change of the particle distribution as the difference be- 
tween the incoming and outgoing fluxes at location r. 

As discussed in the introductory section, in more complex situations, when the linear 
response hypothesis is no longer valid, one observes non-Gaussian behavior, i.e. the long- 
time dynamics is different from that described by the classical Fokker-Planck equation (or 
the usual advection-diffusion equation). At the level of the mean field description, the 
breakdown of linear response means that the particle motion depends on the occupation 
probability in a non-trivial way. The jump probability then becomes a functional of the 
particle distribution function, and Einstein's equation describing the space-time evolution 
of the particle motion must be generalized in order to account for the functional dependence. 

^When = 1, Eq. (1) is the usual diffusion equation with a Gaussian solution [2] and classical scaling 
(7 = 1). 

^In Einstein's formulation the particle jumps arc restricted to symmetrical displacements, i.e. P_|_j = P—j. 
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So introducing Pj {£; t) = pj F (/ {£6r; t)) , with j ^ and where pj is a given distribution 
of displacements (e.g. Pj oc cxp — in (4), we obtain the generalized master equation 

f ir;t + St) - f {r;t) = [^^ (/(r - j<5r; t)) / (r - j<5r; - F (fir^t)) f ir;t)] . (5) 

Generalized diffusion equation. — Along the same lines as the classical diffusion 
equation is obtained from Einstein's master equation, we perform a multi-scale expansion 
of the generalized master equation (5) using an expansion of the time and space derivatives 
of the form 



d_ 

dr 



dr dr 



and a corresponding expansion of the distribution as 
To first order, we have 



0{e') : 
and to second order 



dt 



Sr\ 5(1' 

/o (r; t) = - ( Ji- j — F (/o (r; t)) fo (r; t) 



_ I T^-L\^( "^aF (g) 

^ UtJ dr [ dg 



a(i) a(2) 1 5(1)2 



g=f{r;t) 



/i(r;t)|-( Jji] ^F{fo{r;t))fo{r;t) 



5t 



J2 



(6) 
(7) 

(8) 



-Fih{r;t))fo{r;t) , (9) 



where J„ denotes the moments J„ = ^j=^o j"'Pj- Resummation of these results (see [13] for 
details) yields the hydrodynamic limit of the generalized master equation 



p{r;t)+C-^[Fif{r;t))f{r;t)] = 
D^{Fifir;t))f{r;t)]+l{C'St)-^^E{r;t). 



(10) 



This result is the generalized diffusion equation where C and D are the drift velocity and 
the diffusion coefficient respectively 



St 



(11) 



and 



E(r;0 = _[F(/(r;i))/(r;t)] 



dgF (g) 



dg 



g=f(r;t) 



dr 



[F{f {r-t))J {r-t)] 



(12) 



When there is an external force acting on the particle, Eq.(lO) is further generalized as 
discussed elsewhere [13]. 
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Scaling solution. — Under which conditions is there a scaling solution to the gener- 
alized equation? For simplicity consider the diffusion equation with no drift, i.e. the jump 
probability is space symmetrical and consequently the first moment Ji — and C = 0, and 
(10) reduces to 

§-^fir;t)^D-^[F{fir;t))f{r;t)]. (13) 

Assuming that / (r; t) ~ t (f/t'^) — {^): ^-^id expressing the time and space deriva- 
tives in terms of x, Eq.(13) is rewritten as (see [13] for details) 



-7^a;0 (x) = Dt^-'^^F {r^c/, (x)) 4> (x) 



dx^ 



(14) 



The time-dependence on the right can only be eliminated if F{g) = g'' = t (j/i for some 
number ry, and hence 1 = t^^'^'^t^^^ , i.e. 



7 



1] 



(15) 



Thus, when 7^ 0, this describes anomolous diffusion: (r^) ~ t^+'^ . In this case, Eq.(14) 
becomes 

and admits a q-exponential solution (see [13] for details) 



1 + 2t] 

1+7? 



B 



1 



V 



2 (2 -h 7?) (1 + 27?) BD^ 



(17) 



where B is an integration constant. With rj = 1 — q, and returning to the original space and 
time variables, (17) takes the canonical g-exponential form 



fir;t) 



Bqt — 



1 - (1 - q)Mg 



with 



B„ 



-2(3-<7) {3-2q)B. 



(18) 



(19) 



So, with no drift and no external field, the generalized random walk model describes anoma- 
lous diffusion with q-distributions It also follows from (18) that the distribution of the 
values taken by /(r; t) at any fixed value of time has the form of a power law [14]; 



drD—-t-—5{f{r;t)-f) 



f- 



1 - 



(20) 



An important result of the present analysis is that the power law dependence of the 
transition probability, Pj = pj F{f) with F{f) = is not introduced as an ansatz, but 
follows from the demand for a scaling (or self-similar) solution to the generalized diffusion 
equation. 

Now introducing the power law dependence F{f) — /'' (with 7y > for normalization 
Pj = 1) in the generalized equation (10), we obtain (with 77 = a — 1) 



|/(.;t) + c|:r(.;t) 



Qj.2 • 



(21) 



^One verifies straightforwardly that for q ^ 1, one retrieves the classical Gaussian distribution. 
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Fig. 1: (Color online) Generalized diffusion with no drift: the distribution function /(r; t = 2000 
time steps) obtained from Monte-Carlo simulations (symbols) and the solution of the generalized 
diffusion equation (18) (solid lines) for a = 3 (upper curve) and a = 1 (10''' walkers; initial condition: 
f{r,t = 0)=S{r)). 



with 

E(r;i) = (l-ar-i(r;t))|-r(r;i) . (22) 

Comparison of Eq.(21) and Eq.(l) shows that the two equations are the same in the absence 
of drift (C = 0). Even with nonzero drift, this corresponds to a generalized porous media 
equation when the second term on the r.h.s. of (21) vanishes, i.e. for St 0. So the 
phenomenological generalized porous media equation is an approximation which can be 
obtained in the hydrodynamic limit from the generalized master equation with a power 
law dependence for the transition probability (and in the absence of external force [13]). 
However Eq.(21) contains an additional term which, in general, cannot be neglected (see 
next section). 

Microscopic simulations. — Monte-Carlo simulations are performed with the gener- 
alized master equation (5) using the power law dependent jump probabilities and prescribed 
Pj distributions: pj ~ ^ for j = [—2, +2] (space symmetrical jumps and, so, C=0) and 
Pj = for j ~ [— 2,-f2] (space asymmetrical jumps and so with non-zero drift velocity), 
and the results are compared with the numerical solution of the generalized diffusion equa- 
tion (21). Figure 1 illustrates the case without drift for 77 = 2 (a = 3 and q = —1) showing 
perfect agreement between the Monte-Carlo data and the g-cxponential solution (18); for 
comparison the classical Gaussian result (77 — 0, q — 1) is also shown. 

Two examples with drift are given in Figs. 2 and 3 for a = 1.1 {q = 0.9) and a ~ 2 
(g = —1) respectively showing excellent agreement between the simulation data and the 
solution of the nonlinear equation. We also computed the solution of the generalized diffusion 
equation without the extra term E{r;t) for the value a = 2; the results are given by the 
dashed lines in Fig. 3. The systematic discrepancy with the simulation results gives clear 
evidence that the term given by (22) in the generalized equation (21) cannot be neglected. 
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Fig. 2: (Color online) Generalized diffusion with drift: the distribution function /(r; t = 100 to 
1600 time steps) obtained from Monte-Carlo simulations (symbols) and the numerical solution of 
the generalized diffusion equation (21) (solid lines) for a = 1.1, i.e. q — 0.9 (10''' walkers; initial 
condition: f{r,t = 0) = S{r)). 

To the best of our knowledge the present results provide the first microscopically based 
demonstration of the nonhnear diffusion equation. Further results, including the case where 
the transition probabilities have full spatial dependence (i.e. not only on the distribution 
at the originating location) and the generalization with an external force (i.e. the nonlinear 
advection-diffusion equation) are discussed in [13]. 

* * * 
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Fig. 3: (Color online) Same as Fig. 2 for a = 2, i.e. q = —1. Dashed lines: see text 
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